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^ Abstract: The moment of inertia of a trapped superfluid gas of atomic Fermions 

^ ( Li) is calculated as a function of the temperature. At zero temperature the moment 
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O of inertia takes on the irrotational flow value. Only for T very close to Tc rigid rotation 
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is attained. It is proposed that future measurements of the rotational energy will 



unambiguously reveal whether the system is in a superfluid state or not. 
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1. Introduction 

The advent in 1995 of Bose-Einstein-Condensation (BEC) of atomic Bosons in 
magnetic traps certainly represents a milestone in the study of bosonic many 
body quantum systems. This is so because a systematic study of these systems, 
starting with the free particle case, as a function of increasing density, particle 
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number, and other system parameters seems possible and has already 
progressed to a large extent while going on at a rapid pace [1 ,2]. On the other 
hand the recent experimental achievement of trapping ®Li atoms and other 
fermionic alkali atoms [3] also spurs the hope that for the fermionic many body 
problem, as much progress will be made in the near future as for the bosonic 
systems. Indeed the first Fermi-Dirac degeneracy of trapped '^°K atoms has 
already been observed (see B. De Marco and D.S. Jin [4]). In this reference also 
more of the physics of trapped Fermionic atoms is discussed. For atoms with 
attractive interaction one can envisage that the trapped system undergoes a 
transition to the superfluid state. For instance ^Li atoms can be trapped in two 
different hyperfine states. In the spin polarized case the s-wave interaction turns 
out to be very strong and attractive (scattering length: a = -2063 ao with ao the 
Bohr radius) favoring a phase transition to the superfluid state. This possibility has 
recently provoked a number of theoretical investigations (see [5] for a more 
detailed discussion of a possible superfluid state). One major question which is 
under debate is how to detect the superfluidity of such a fermionic system, since 
in contrast to a bosonic system the density of a Fermionic system is hardly 
affected by the transition to the superfluid state [6]. Several proposals such as the 
study of the decay rate of the gas or of the scattering of atoms off the gas have 
been advanced [5]. Though such investigations may give precious indications of a 
possible superfluid phase, we think that in analogy with nuclear physics, a 
measurement of the moment of inertia certainly would establish an unambiguous 
signature of superfluidity. To measure the spin and the rotational energy of 
trapped atoms definitely is a great challenge for the future. However, in nuclear 
physics, where y-spectroscopy is extremely well developed, the strong reduction 
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of the moment of inertia witli respect to its rigid body value lias been considered 
as a firm indicator of nucleon superfluidity immediately after the discovery of 
nuclear rotational states almost half a century back [6]. Therefore awaiting future 
experimental achievement also for trapped fermionic atoms, it is our intention in 
this work to give some theoretical estimates of the moment of inertia as a function 
of deformation of the traps or temperature of the gas. In this study we can largely 
profit from the experience nuclear physicists have accumulated over the last 
decades in describing such phenomena. The expectation is indeed that there will 
be a great analogy between the physics of confined atomic Fermions and what 
one calls in nuclear physics the liquid drop part of the nucleus. As astonishing as 
it may seem assemblies of fermions containing no more than ~200 particles 
(nucleons) already exhibit an underlying macroscopic structure well known from 
the Bethe-Weizsaecker formula for nuclear masses [6]. In superfluid rotating 
nuclei as early as 1959 Migdal proposed a statistical description of the nuclear 
moment of inertia [7] which grasped the essential physics of a self contained 
rotating superfluid Fermi liquid drop and which serves as a reference even today. 

In the present work we will cast Migdal's approach into the more systematic 
language of the Thomas Fermi theory which together with its extensions is 
applied since long to normal fluid but also to superfluid nuclei [6,8,9]. It is 
fortunate that we can profit from this experience for the description of trapped 
fermions, since their number of order 10^, together with the smoothness of the 
potential, certainly turns a statistical description into a very precise tool. On the 
other hand it may not be excluded that in the future the study of much smaller 
systems of trapped atomic Fermions with numbers -10^ may be studied probably 
revealing many analogies with nuclei such as shell structure etc. The investigation 
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of the transition from microscopic to macroscopic as tine number of particles is 
increased continuously may then become a very interesting field also in the case 
of atomic Fermions. In detail our paper is organized as follows. In sections 2 and 
3 we review the Thomas Fermi approach to inhomogeneous superfluid Fermi 
systems. In section 4 first the so-called Inglis part of the moment of inertia of a 
rotating superfluid and confined gas of atomic Fermions is presented. Second the 
influence of the reaction of the pair field on the moment of inertia is calculated. It 
is shown that this leads to the irrotational flow value in the limit of strong pairing. 
In section 5 the current distributions in the superfluid and normal fluid regimes are 
contrasted. In section 6 the numerical results are presented in detail together with 
discussions and conclusions. 



2. Thomas Fermi approach to fermionic atoms in deformed traps 

The Thomas-Fermi approach (TF) to trapped gases of atomic Bosons is a well 
accepted practice [2]. For trapped atomic Fermions the same approximation 
underlies different conditions. It has, however, recently also been applied to this 
kind of situation [5]. The TF approach for Fermions also is extensively applied to 
other finite systems such as atomic nuclei, metallic clusters, etc. The smallness 
parameter is given by 



Where V is the mean field potential and kp{r) the local Fermi momentum. 
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(2.2) 



With a typical frequency of tine external harmonic potential of cOq=7 nK and 
Ep = 600 nK one realizes that ;/ « 1 up to close to the end of the classically 
allowed region. For integrated quantities the region around the classical turning 
point carries little weight and therefore the TF approximation for a number of 
trapped atoms of the order of 1 0^ is certainly very well justified. 

Furthermore, as in the boson case the TF approach [6] to trapped atomic 
gases becomes extremely simple by the fact that the large interparticle distance 
makes a pseudopotential approximation to atomic interactions valid. Let us 
therefore write down the TF equation for a doubly spin polarized system of 
trapped (^Li) atoms in the normal fluid state. For convenience we first consider 
the system at zero temperature T discussing the T ^^^^ case later on. In TF 
approximation the distribution function for particles in each spin state is given by 
(in this work we only consider equal occupation of both spin states) 



f(R,p) = 0{M-Hj, 



(2.3) 



with 



(2.4) 
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Here n is the chemical potential, and v;,;,(R) stands for the trap potential which is 
supposed to be of harmonic form. The density p(R) is obtained from the self- 
consistency equation, 




(2.5) 



with. 



^,(R) = ^|^U-yjR)+gp(R)), 



(2.6) 



the local Fermi momentum. The coupling constant g is related to the scattering 
length in the same way as in the case of Bose condensed gases [1 ,2 ] via 



The TF equation (2.5) leads to a cubic equation for the self-consistent density, 
which straightforwardly can be solved as a function of the external potential. In 
this paper our main interest will be the study of the moment of inertia of a rotating 
condensate. Since the study is very much simplified assuming that the self 
consistent potential is again a harmonic oscillator and since the effect of the 
attractive interaction between the atoms essentially results in a narrowing of the 
self consistent potential with respect to the external one we will use instead of the 
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exact TF solution for the density tine following trial ansatz for the local Fermi 
momentum 



(2.8) 



where <w^,<Wyand <w^are the variational parameters. The chemical potential is 
determined from the particle number condition 



(2.9) 



and the kinetic energy density is given by 



(2.10) 



We then can analytically calculate the total energy 



(2.11) 



as a function of <w^,<Wyand co^. Minimizing this expression with respect to 
0)^,0) and cojor a given external deformed harmonic oscillator potential, 



(2.12) 
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leads to the variational solution. For the spherical case (o^ = (o^ = (o^ = (o, this is 

shown in Fig. 1. We see that this approximation to the TF equation is quite 
reasonable. For an external harmonic potential with frequency v = 6)o/2;r = 144Hz 

or fKo^jkj^ =6.9nK, corresponding to the conditions of the experiment of Bradley 
et al. [10], the variational frequency is %(ojkj^ =7.69nK. Since (oxo^ this implies 
a compression of the density. Increasing (o by 6% (^ftj/A:^ =8.21nK) from its 
variational value allows an almost perfect reproduction of the full TF solution. We 
will adopt this latter value in all our forthcoming calculations. The experimental 
situation for the rotating deformed case is such that the rotation of the trap is 
performed around the x-axis (the long axis), permitting slight triaxial deformations 
in the plane perpendicular to the x-axis that is in the y-z plane. In order to simulate 
such an experimental situation we simply first make a volume conserving 

{(o^)'(o^ = (o\ (2.13) 



prolate deformation around the x-axis, 

o = ^ (2.14) 

CO, =ft)<7~'^^ ; 0). = 0)= Q)^ 

In order to increase the central density there is interest to make rather strong 

eccentricities and cojcoj^ =- is a typical value which we will adopt in this paper. 
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In a second step we fix co^ and deform around the x-axis again keeping tine 
volume fixed. We define the deformation parameter as 



5 = ^ 



(2.16) 



We finally have the two-parameter deformation 



2/ 

0)^ = ^ycj"*^^^"*^' 



(2.17) 



with <7 « ^ < 1 . From now on we therefore will use for the nonsuperfluid Wigner 
function at zero temperature the expression 



/(R,P) = ^ 



2m 



(2.18) 



with 



(2.19) 



and Q}^, cOy, co^ from (2.17) and // determined from the particle number 



condition. 



3. The superfluid case 
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Since, as described in tine introduction, trapped spin polarized ^Li atoms, in 
different hyperfine states, feel a strong attractive s-wave interaction, the system 
very likely will undergo a transition to the superfluid state at some critical 
temperature Tc as was discussed in detail in ref. [5]. As we have pointed out in 
the introduction the superfluid state will unambiguously reveal itself in its value of 
the moment of inertia. At the moment the measurement of angular momenta of 
trapped Bose or Fermi gases has not been achieved and represents a future 
challenge to the experimenters. In order to establish how the two essential 
system parameters which are the value of the gap, the temperature of the 
system, and the deformation of the external trap influence the value of the 
moment of inertia, we will now proceed to its evaluation in the superfluid state. 

Since we are dealing with an inhomogeneous system, even in the 
nonrotating case the gap is actually a nonlocal quantity A(r,rj or in Wigner 
space A(r,p). We will find later that at zero temperature the coherence length of 

fi^ k 

the Cooper pair, ^ = ^, is larger than the oscillator length 



careful with applying the TF theory for temperatures T much lower than the 
critical temperature Tc where the gap vanishes. We will discuss this point more 
thoroughly later and in the appendix. We therefore go on and apply the TF 
approximation to the superfluid state. It has been shown in [9,6] that to lowest 
order in h the gap equation is given by 



m A 



1 = ^. 



'meoQ ~ 0.63x10^ flg with the Bohr radius. We therefore have to be 




A(R,k) 
2£(R,k) 



tanh 



V 




(3.1) 
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where ^'(R,?) is the quasiparticle energy, 



2 

Pf 



(R)T 



2m* (r) 



+ A^(R,p) 



(3.2) 



with pp{R) = hkp{R) the local Fermi momentum (2.6). Since the effective mass 

m* is so far unknown for trapped gases of atomic Fermions we will take m* = m. 
Furthermore, for the time being, as in [5], we will eliminate the interatomic 
potential v in (3.1), expressing it by the scattering length (2.7). We then obtain [5] 



tanh 



E{R,k) 



2T 



2£(R,k) 2{e,-£p{R)) 



A(R,k), 



(3.3) 



where P stands for principal value, = — — and = — ^ . At zero 



2m 



2m 



temperature, as described in [11], (3.3) can be solved analytically in the limit 

A(R,p^(R)) 



^ . The result is given by 



A, (R) - A(R, kp (R)) = Se-'ep (r) e ^^^^^^'H 



(3.4) 



A posteriori one can verify that — « 1 for all values of R and therefore (3.4) is 
an excellent approximation to (3.3). This also has been found in [5]. For 286500 
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®Li atoms, the case considered in [ 5], tine gap is sliown for a splierical trap as a 
function of the radius in Fig. 2. 



For the determination of the critical temperature Tc and, later on, for the moment 
of inertia we will need the value A of the gap at the Fermi energy. Since the 
detailed level structure at the Fermi energy is unknown and in fact unimportant, 
we will consider the gap A{sp) averaged over the states at the Fermi energy 

A{£,)^A = Tr(Ap{e,)), (3.5) 



with 



Pis) = 4l 1 1 1^^"^ -£j = ^M£-H), (3.6) 

where |n)and £„ are the states and energies of the harmonic oscillator with 
frequency co and 

g{£) = Y^d{6-£j = Tr5{s - H) , (3.7) 

n 

is the level density. 

It has been shown in [12] that again the TF approximation leads to an 
excellent average value 

g {£p)^ [iTm] 
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In the spherical case with A^(i?) from (3.4) all integrals but the radial one can be 
performed analytically, the latter being done numerically. For the case shown in 
Fig. 2 one obtains 



A = 16.4 nK. 



(3.9) 



Quantum mechanically the BCS equations should be solved in the self-consistent 
Hartree-Fock (HF) basis and then Tc is a global parameter which must be 
determined from the solution of the quantum mechanical gap equation. Since we 
believe that the value in (3.9) comes rather close to the quantum mechanical 
value of the gap at the Fermi energy, we can obtain Tc from the usual BCS weak 
coupling relation [13] A = 1.76r^ to be 

r^^lOnK. (3.10) 



nK which corresponds to our approximate "self-consistent" harmonic solution with 
fy = 8.21nK and kp\a\ = 0.56 one obtains ^ = 4.10^ which is about a factor of 

seven larger than the oscillator length of the trap (see above) which contains 
286500 particles. This seems to invalidate the TF-approximation. We, however, 
know by experience that often the TF-approximation remains quite reasonable 
beyond its limit [6]. For example the conditions of validity in [9,12] for superfluid 
nuclei are much worse than here and still the results are accurate beyond 



From (3.9) we obtain the coherence length ^ 



ft kp ^ 

m A ^ kp 



. With £p = 983.67 
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expectation. We therefore think that the values (3.9, 3.10) are reasonable 
estimates for the gap and the critical temperature. In order to check this 
assumption we give in the appendix a more refined semi-classical solution of the 
gap equation which only demands that the TF-approximation in the normal fluid 
state is well justified. We find values for A and Tc which are ~ 30% lower than in 
(3.9)(3.10). In view of the crudeness of the TF-approach this indicates a quite 
satisfying consistency between the results. 

We also will have to know the detailed T-dependence of the gap A(r) which 
however, in BCS theory, given A(o) and Tc , is determined by the universal 

A(r) T 

function — W^in terms of — . This function is determined from the solution of the 

A(o) t; 

equation [13] 



■In 



A(0) 



= A 



T 



with 



(3.11) 



A{u)=\dy^= 



f f 1—; 



l-tanh 



)) 



For completeness it is shown in fig. 3. This T-dependence of the gap we will later 
use for the evaluation of the moment of inertia. 



4. Moment of inertia 
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The moment of inertia of a rotating nucleus has fully been formulated in linear 
response theory (i.e. RPA) by Thouless and Valatin [14]. The corresponding 
expression is therefore called, in the nuclear physics literature [6], the Thouless- 
Valatin moment of inertia. It consists of two parts, the so-called Inglis term, which 
describes the free gas response, and the part, which accounts for the reaction of 
the mean field and pair potential to the rotation. In the superfluid case the Inglis 
part has been generalized by Belyaev [15] and the linear reaction of the gap 
parameter onto the value of the moment of inertia was first evaluated, together 
with the Inglis term, by Migdal [7]. The reaction of the HF field on the rotation is a 
minor effect and we will neglect it in this work. We therefore will write the moment 
of inertia as a sum of the Inglis-Belyaev term 0^_g and the Migdal term 0^ . In 
total 

= 0,_^ + 0^. (4.1) 

In order to derive an expression for in linear response theory we will use 
the Gorkov approach described in detail in many text books (in what follows 
we will use the notation of [16]). Since in addition the derivation of the linear 
response for is given in the original article of Migdal [7] and rerepresented 
in a more elaborate version in [8], we will be very short here and only give 
more details where in our opinion the presentations in [7,8 ] may not be 
entirely explicit. Let us start writing down the Gorkov equations in matrix 
notation 



15 



f 3 

- — -H + ju 



G = l-AF' 



(4.2) 



+ H* +JU 



= A G 



with 



H — Hq — = H Q + , 



(4.3) 



where now is the shell model Hamiltonian (2.4 ) or rather the approximate 
one used in (2.18), (2.19) and 

L =r p —rp , 

the angular momentum operator corresponding to a rotation with angular 
frequency £i around the x-axis. In (4.2) G and F are the normal and anomal 
Matsubara Green's functions (see Chap. 51 of [ 15]) 



(4.4) 



F%„'=-(r,a/(r)a/(r')' 
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Linearising (4.2) witli respect to H^, tliat is G = Gq + G^, F"^ = + and 
A = Aq + Ai (as mentioned we will neglect the influence of the rotational field on 
Hq) one obtains for (4.2) 

G.=Gu-B+Gm, (4.5) 
with 



GiM - -Gq^Fq'' - Fq'^Ai Go 



(4.6) 



and 



F; = D,H;F,' + F.'Hfi, - F,'A,F,' + D^'Go ■ (4-7) 



where, 



and are the Matsubara frequencies [16]. 

In (4.5,4.7) we have split the first order Green's function in an obvious notation 
into the Inglis-Belyaev and Migdal contributions. For the latter one needs the 
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linear reaction of tine pair field to the rotation. We later will see how this can be 
determined from (4.7). First let us, however, evaluate the l-B part of the 
moment of inertia. 



4.1 The Inglis-Belyaev part of the moment of inertia 

The l-B part of the moment of inertia can be evaluated without the 
knowledge of i.e. without the use of (4.7). The density response 

corresponding to Gj^.^ of (4.5) is evaluated from the limit t'^t^ or from 
summing over the Matsubara frequencies in the upper half plane (see Ch. 7 of 
[16]) . One obtains the well known result [6, 7, 8, 15,16] 

(A/-bL = («|^>X'' (4-8) 



with 



F = F,(l-/-/j+F(/-/j, (4.9) 



where 



F,=F.fe..j= ^-^^;;y;.-^j^:'f--> , (4.10) 

n n\ n — n f 



/ = /k)=^^ ; /'=/k') (4.11) 
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and 



(4.12) 



are the quasiparticle energies with as ingredient s„, the energies of the 
harmonic oscillator potential (2.19). The gap parameters A^have been 
replaced in (4.10), in analogy to (3.5), to statistical accuracy by a(^J, the 
ones averaged over the energy shell. The moment of inertia is given by 



Since we are interested at temperatures T<T^, which are very low with 

respect to the Fermi energy, we checked that one can to very good accuracy 
neglect in (4.9) the thermal factors (4.11). The only important temperature 
dependence of the moment of inertia therefore exists via the T-dependence of 
the gap. We thus will henceforth treat all formulas as in the T=0 limit keeping, 
however, the T-dependence of the gap. With this in mind we can write for the 
moment of inertia 



In this formula the important quantity to calculate to statistical accuracy is, 



(4.13) 







I-B 



(4.14) 
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L/ (n, n') ^ \(n |L, | n')|' = rr[(L, ) ( | n '> (n'|L, | n)(n \ )] 



(4.15) 



where = o(R,p) means the Wigner transform of the operator O [6]. To this 
purpose we again replace the density matrices \n)(n\ and \n)[n\ by their 
average on the energy shell (3.6) 

\n){n\^p{e„). 

We therefore obtain 



with 



(4.16) 



{lXco,coX^[5[(o'-H,]l^8(co-H,)\ 



(4.17) 



Introducing into (4.17) the Fourier representations of the two ^-functions and 
transforming to center of mass and relative coordinates one obtains 



with 



J. CO + CO' 

E = ; £ = (0-0) 



and 









-iH„T J 




-iHr.T 

e " 


.2. 





(4.18) 



w 



(4.19) 
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O{t) = e'"''O{0)e-'"''' . 



(4.20) 



To lowest order in ^iwe replace the triple operator product in (4.18) by the 
product of their Wigner transforms [6] 

















-iH„T 








.2. 






.2. 



lim 



and therefore 



(L, {co, coX =L^{E,e,R,p) = U{E- H,^, )\ L 



1 



v2y 



(4.21) 



(4.22) 



with 

Lt{t) = R^{t)p^{t)-R^{t)p^{t). (4.23) 

At this point the choice of our approximate self consistent potential of 
harmonic oscillator form (see 2.19) turns out to be very helpful, since the 
classical trajectories in (4.22) can be given analytically 



R^{t) = cos{h(Oit)+^^cos{hQ).t) 
mcOi 



Pi {t) = Pi cos{h(Oft) + mcOjRf cos{h(Ojt) 



(4.24) 



with i = x,y,z ■ 

In the phase space integral of (4.15), for reasons of symmetry, only the 
diagonal terms of L^'.L^'(^) survive and therefore we obtain 
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+ 



Q) 



2^y . „ 2 <W, 



0), 



CO,, 



smi^co t)sm{hcojt) 



where 



1 



2m 



-,3/ 



{E-V) 



(4.25) 



(4.26) 



is the density in TF approximation (see eq. (2.5)). 

The product of cosine and sine in (4.25) can be expressed in terms of the 
cosine of the sum and difference of the arguments and then the t -integral in 
(4.23) can be performed. This leads to ^-functions which allows to perform 
also the ^ -integral. Furthermore, as shown by Migdal [7] 



^ e ^ 



F,[E,£)^ 1-G 
where (see eq. (3.5)) 



5{E-ii), 



(4.27) 



^=^{ep)., 



and 



G{x): 



arsi 



nh(x) 



(4.28) 



Finally one obtains for the l-B part of the moment of inertia the following 
analytical expression [ 7,8] 
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=0 

^I-B ^ri. 



rigid 



2 2 

0) +0)^ 



(4.29) 



where 



(0+ = cOy±o)^ ; G+ = G 



and 



(4.30) 



(4.31) 



is the moment of inertia of rigid rotation. From (4.29) we see that 



lim0 =0 ■ lim0 = 

\ n I-B rigid > \^^^^^I-B ^ 

A— >0 A— >oo 



The latter result is clearly unphysical and we will see how the account of the 
reaction of the pair field on the rotation will reestablish the physical situation. 

4.2 The Migdal term 



The density response corresponding to the Migdal term is obtained from 
(4.6) 



VMmW 2eeIe +e.) 

n n \ n n / 



(4.32) 
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In (4.32) we need to know which we can gain from (4.7) in the following 
way; in the limit z' we obtain from F/the anomal density , 



\ i /nn 



IE E ,(E +E .) 



(4.33) 



In analogy with the non-rotating case where /Tq = — > we also have 

2,E 



iic:) . = -A* 

\ 1 /nn In 



1 1 

- + - 



4E 4E' , , 



(4.34) 



This relation stems from the fact that the quasiparticle energies contain the 
gap only in the form AA* and therefore there is no further first order correction, 
since in our case the external field is a time odd operator and thus 



(4.35) 



Equating (4.33) and (4.34) yields 



KKn'\n' + 2Ao„//,„„.^„, + 2A„„A,„„,A„„, + (aI„ + A^, + (^„ - ) A^^ 

2E„E„,{E„+E„,) 



= 0. 



(4.36) 



At this point we again exploit the fact that expression (4.36) is strongly peaked 
around the Fermi energy surface. Following [7], in analogy with (4.27), we have 



e„ -e 



^ 5 



V 2A , 



n n 



(4.37) 
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With (4.35) we then obtain for (4.36) 



{n\L\n\ 



2A 



+ 



2A 



X n 



£-£ 



2A 



£ +£ , 

n n 



= 0, 



(4.38) 



where Stands for the time derivative ofL^. Summing on n and n' and following 

exactly the same line of semi-classical approximations as the ones used for the 
derivation of one arrives at the following relation [8] 



jdtG{T)\ 



p 

[iTthf 



2A 4A' 



(4.39) 



where g(t) is the Fourier transform of g(x) (4.28). 
For the potential in (2.19) (4.39) is solved by 

x{^) = aR^R,, (4.40) 
with 

G + G 

a = -2AmO)^(0 ^ — . (4.41) 

^ ~(olG^ + olG_ 

Inserting this solution into (4.32) leads for the Migdal part of the moment of inertia to 
[7,8] 

- ^ rigid ^9 , ^,2 ^,2^ , ^,2^ ■ \^-^^) 

CO^ + C0_ CO^G^ + (0_G_ 



Together with (4.29) the expression for the moment of inertia is now complete. Let us 
again mention that we neglected the temperature dependence besides the one 
contained in A = A(r), since all other T-dependence for T<T^ is negligible. The 

moment of inertia can then be calculated as a function of deformation and 
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temperature. For example it is immediately verified that for A ^ oo (4.42) yields the 
irrotational flow value 




(4.43) 



and therefore 



lim © = lim + 0^ ) = © 



irrot ' 



(4.44) 



which is the correct physical result. 
5. Current distribution 

Other quantities, which may be interesting also from the experimental point of 
view, are the current distributions of the superfluid rotating gas. Indeed after a 
sudden switch off of the (rotating) trap the atoms will expand keeping memory of their 
rotational state. So if the velocity distribution of the expanding atoms can be 
measured, one may be able to deduce the rotational motion the atoms have had 
before the trap was taken away. The current distribution, as we will see, depends, as 
the moment of inertia, strongly on the superfluid state of the gas. In order to calculate 
the current distribution we first write down the Wigner function of the density 
response which can easily be read off from the formulas given in Section 4. In 
obvious notations we obtain [8] 
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-Q. 



(O^G_ - 0)_G^ (O^G_+(0_G^ 

-RyPz -RzP 



2a), 



y 



(5.1) 



a 



Am(R.p) = — 

2m A 



co.G. - co G „ 6), G, +0) G „ 



Siju-H,J, (5.2) 



With the usual definition of the current 



one obtains 



(5.3) 



^ co^G_ + (0_G^ 
2a)„ 



(5.4a) 



^_(O^G_-(0_G^ 
2(0, 



(5.4b) 



(0+(0_ (G_ + ) (0^G_ + co_G^ 
coIg^ + co^G CO,, 



(5.5a) 



co^co_ (G_ + G^ ) (o^G_ - co_G^ 
(olG^ + (olG_ (O^ 



(5.5b) 



with of course, =0. 
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Again we see that in tine limit A ^ oo tine current approaclies tine correct irrotational 
flow limit 



A->oo OJ 




2 2 



(5.6) 



whereas in the limit of A ^ we obtain a rigid body current. As we have seen for , 
as a function of temperature and deformation, we easily can go from one limit to the 
other. 

6. Results and conclusion 

We show in Fig. 4 a and b the current distribution for the two extreme cases of 
irrotational and rigid body flow in the laboratory frame respectively. We see that the 
flow pattern is completely different in the two cases. In Fig. 4b the flow pattern clearly 
corresponds to rigid rotation of an ellipsoid with the long axis in z-direction. Also Fig. 
4a represents a typical irrotational flow pattern well know from hydrodynamics. As a 
function of temperature one continuously can pass from one flow pattern to the other. 
The point we want to make is that for small deformations S , as can be seen from 
(5.6) there is almost no irrotational current for low temperatures and this will then be 
reflected in a very low rotational energy as we will discuss now. 

In Fig. 5 we show as a function of a{t) and with Fig. 3 also as a function of 

T. We see that for a typical eccentricity ^ = ^ = 0.8 the moment of inertia changes, 

as a function of temperature, by large factors. At r = the gap values found in this 
paper are in the range of 10-20 nK and therefore the moment of inertia is close to its 
irrotational flow limit. This actually means that the moment of inertia is, with respect to 
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its rigid body value, very small, since for ^ ^ 1 , i.e. for spherical symmetry around 
the rotational axis (x-axis) the moment of inertia goes to zero (see (4.43, 4.44)). 
Consequently in this case the gas is not at all following the rotation of the trap. 
However, increasing the temperature i.e. decreasing the gap value has a dramatic 
influence on the rotational motion of the gas, since in the range < T < Tc the 
moment of inertia raises very steeply attaining its rigid body value for T = Tc. In this 
limit the gas rotates as a whole with the same angular frequency as the trap. The 
abruptness of the raise is the more pronounced, the smaller the eccentricity S is. 
(see Fig. 5). Experimentally non-destructive or expansion imaging can be used to 
watch the gas rotate and then the rotational energy 

Ero.=^^\ (6.1) 

can be obtained by integrating the angular velocity over the density profile *^ The 
rotational energy therefore directly follows the variation of the moment of inertia. One 
deduces that the measurement of the variation as a function of T of the rotational 
energy should be well within the experimental possibilities, once the technique of 
putting the trap into rotation has been mastered. 

In our discussion we have ignored the possibilities of vortex formation. The 
determination of the onset of instabilities versus vortex formation in a finite Fermi 
system is not a completely easy task and we will postpone such an investigation to 
future work. 

* ) We are grateful to the referee for pointing to this possibility. 
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However, since the rotational frequencies Q. considered in tliis paper are mucin 
smaller than the oscillator constant coo (il/coo « 1) we think that our result will not be 
spoiled by the appearance of vortices. An indication can also come from the case of 
trapped Bosons where vortices, depending somewhat on the number of atoms, do 
not appear for values n/coo < .5 (see ref. [2]) 

From the above discussion we see that it may well be in experimental reach to 
reveal an eventual superfluid state of the gas once the technique of putting the trap 
into rotation will be mastered experimentally. A closely related phenomenon to 
rotation is the so-called scissors mode which was originally discussed and found in 
deformed nuclei [17] and then proposed [18] and also very recently found [19] for 
trapped Boson condensates. Suppose the trapped atomic system is rotating very 
slowly and suddenly the rotation of the deformed trap potential will be stopped. Due 
to inertia the atomic cloud will continue rotating back and forth around the fixed trap 
position if the initial rotation was gentle enough. If for the purpose of a rough 
argument we suppose that this oscillatory motion has so small amplitude that in a first 
approximation we can neglect shape distortions of the cloud, then, if the oscillations 
are in harmonic regime, the frequency of the scissors mode is given by 



where C is the constant of the restoring force. The frequency 0)^ will strongly depend 

on whether the system is in the superfluid state or not. In this way the above cited 
experiment has indeed unambiguously revealed that the Bose condensate is in a 




(6.2) 
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superfluid state [19]. It is evident tliat scissors modes could also be excited in trapped 
Fermi systems, as this was already mentioned in [18]. Since in Fermi systems for 
temperatures T~Tc one can suppose that the temperature dependence of the force 
constant is weak with respect to the one of the moment of inertia 0, one will find a 
strong difference for the value of 0)^ in the superfluid and unpaired regimes 

respectively. As long as the temperature is so small that the normal fluid component 
can be neglected, the temperature dependence of co^ can be deduced from the one 

of obtainable from Figs 3 and 5 of this work. We are, of course, aware that the 
experimental situation may be more complicated needing a more refined discussion 
similar to the one given in [18]. A more detailed investigation of the scissors mode for 
trapped Fermions may be given in future work. 

In summary we proposed in this work to measure the dynamics of a rotating trapped 
gas of atomic Fermions as a function of temperature and deformation to detect 
whether the system is in a superfluid state or not. Quite detailed and quantitative 
calculations for the moment of inertia and velocity distributions have been presented. 
Other quantities well studied in the case of rotating superfluid nuclei [6] such as Yrast 
lines, even-odd effects, particle alignment, etc., may also become of interest in this 
case. 
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Appendix 



In this appendix we want to give a more refined semi-classical solution of the gap 
equation. Let us write the quantal version of (3.3) at T = in BCS approximation [6] 



where n labels the states of the (spherical) harmonic oscillator with single particle 



energies and n is the time reversed state. As usual E^=-^{£^-£pf +A\ is the 

quasiparticle energy and (nn\\\n'n'y\s the matrix element of the zero range two body 
force 



Since what matters is the gap at the Fermi level and since for particle numbers of the 
order 10^ the degeneracy of the oscillator shells is very high it seems a very 
reasonable approximation to replace in (A.1) all quantities by their corresponding 
values averaged over the energy shell (3.5, 3.6). Equation (A.1) can then be written 
as 



tanh 




(A.1) 




(A.2) 



tanh 




(A.3) 
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where gie) is the level density (3.7) and v(f,f') is the averaged two body matrix 
element 



Y{e,e')= I Y^S{e-e„W-£n'){nn\y\n'n). (A.4) 

At this stage one could try to solve the gap equation numerically. However, again in 
view of the huge number of particles it is certainly a good approximation to pass to 
the Thomas-Fermi limit. For the level density this is immediate. The TF limit of 
(A.4) can be obtained in locally summing over plane waves and we obtain 



3 ^ inf(/f,/i) 



where is the classical turning point given by e = v{r^) and v{r) = m(o^r^ jl is the 

harmonic oscillator potential. We have made numerical check that (A.5) is indeed a 
good approximation to the quantal counterpart for the case of large particle numbers 
[20]. We notice that (A.5) only needs the TF approximation in the nonsuperfluid state 
where it is well justified (see sect. 2). Having an expression for average level density 
and matrix element at hand we can proceed to solve (A.3) We will do this again in the 
limit A{£p)/£p «1 and obtain (see [1 1]) at T = 



64k p \a 



+ I[£p 



(A.6) 
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with 



1 v^'' 



\-x 

The integral /(f^) is evaluated numerically and we obtain 



K{e,) = ^e,e-''''e ^'^^^'''1 . (A.8) 



With £j, = 983.67 nK which corresponds to our "self-consistent" harmonic solution and 
yt^ a =.56 one obtains A(£-^) = ii.29nK. This value is about 30% smaller than the one 

extracted in (3.9), which, however, in view of the roughness of the TF approximation 
can be considered as a rather satisfying consistency of the results. 
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Figure captions 



Figure 1 : Density profiles, for tine case of a splierical trap, of tine non-interacting 
case (full line), the interacting case once calculated exactly from (2.5) (crosses) 
with Y^^ given by (2.12) and once using the variationally determined harmonic 

oscillator potential (open squares). Squeezing the variational (o by 6% yields a 
density which lies on top of the exact TF solution. 

Figure 2 : The gap for a spherical trap as a function of the radius 

Figure 3 : Ratio of the energy gap to the gap at T=OK as a fun ction of 
temperature 

Figure 4 a and b : The current distribution for the two extreme cases of 
irrotational (a) and rigid body (b) flow in the laboratory frame. In both cases the 

deformation parameters are set to (<7 = -,^ = .8 ) and the angular frequency £i 

8 

around the x-axis to 1 nK. 



Figure 5 : The moment of inertia as a function of the gap for different values of 

the deformation 5 = — and o- = - . 

(o, 8 
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